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Recent progress in Lattice QCD is liighliglited. After a brief introduction to tlie 
methodology of lattice computations the presentation focuses on three main topics: 
Hadron Spectroscopy, Hadron Structure and Lattice Flavor Physics. In each case 
a summary of recent computations of selected quantities is provided. 

1 Introduction 

Lattice QCD has come of age. What is meant with this statement is that 40 years 
after its inception [l] this framework is now able to deliver clear-cut predictions for 
a number of phenomenologically relevant quantities. 

From an experimentalist's viewpoint what matters is the precision/accuracy of a 
computation; to be relevant it must match the experiment's precision (which often 
is a challenge). From a theoretician's viewpoint the key issue is that all quantities 
are well defined and that the total uncertainty of the final result (usually split into 
statistical and systematic contributions) can be reliably assessed. In that respect 
the last decade has brought a major improvement: The quenched approximation 
(in which the quark loops that come from the functional determinant are omitted) 
is gone. While often found to have a small impact in practice (a posteriori, i.e. 
by comparing to the fuU/unquenched result), this approximation was a nuisance, 
because its impact was extremely hard to quantify. Todays lattice computations 
set new standards, since the collaborations spend an enormous amount of effort to 
control and quantify all sources of systematic uncertainty. 

There is a number of physics questions in which Lattice QCD plays a key role. 
One example is "Do we understand the weight of this world" (the visible matter 
in the Universe) ? With the Higgs particle established at the 4cr level (as reviewed 
at this conference [2|3| ) one might think that this is the cause. Looking at actual 
numbers we realize that the Higgs mechanism is almost irrelevant for the mass of the 
proton; the average of the up and down quark masses is = {m^^^-{m^^^) /2 ~ 

3.5 MeV in (MS,2GeV) conventions, while the proton mass is Mp ~ 940 MeV. 
At this point we stand accused of comparing apples to pears; the former mass is 
a scheme and scale dependent quantity, while the latter one is a truly physical 
quantity. A more meaningful attempt would relate the proton mass to the mass 
that the nucleon would have in the 2-flavor chiral limit, that is to Mjv — 880 MeV 
with rriuci ^ and all other quark masses held fixed (cf. Sec. [4]). Still, the message 
is the same: The weight of the world which we experience is essentially due to the 
relativistic dynamics that is involved in the process by which the strong force binds 
quarks and gluons into protons (and other color neutral objects). 
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To date the lattice is the only known theoretical tool which can "solve QCD" 
in practical terms and with fully controlled systematics. This means that it is 
capable of establishing the link between fundamental parameters of the Standard 
Model (e.g. quark masses or CKM matrix elements) and the quantities measured in 
experiment (e.g. hadron masses or leptonic and semi-lcptonic widths or branching 
fractions). In the following the goal will be to explain the theoretical underpinning 
of these computations, and to highlight some of the latest results. 

The remainder of this review is organized as follows. Section [2] gives an outline 
of the field-theoretic setup of Lattice QCD computations. The next three sections 
contain the core of the presentation, an update on Hadron Spectroscopy (Sec.|3]), 
Hadron Structure (Sec.|4]) and Flavor Physics (Sec. [5], with a brief presentation of 
the recent FLAG activities). To give the reader an idea of what a tiny segment of 
lattice studies has been covered, we shall attempt to list some of the topics omitted 
in Section [6] before presenting a summary in Section [7]. 

2 Lattice Field Theory 

2.1 Regulating QCD in the UV and the IR 

Quantum Chromodynamics (QCD) is a field theory. As such it requires, in an 
intermediate stage, regularization both in the ultraviolet (UV) and in the infrared 
(IR) to be mathematically well behaved. 

The lattice approach (more precisely: the multitude of lattice approaches, see 
below) does this by replacing the spacetime manifold by a regularly spaced 4- 
dimensional grid. Typically, a hypercubic set of grid points x € (71.1,71,2,71.3,71,4)0 is 
used, with Ui (i — 1,2, 3) running from 1 to Ng and 714 running from 1 to Nt and a 
the lattice spacing. The temporal direction is considered the fourth direction, since 
we work in Euclidean space, where all directions have the same sign in the metric. 
The lattice spacing a is usually chosen isotropically, i.e. the same (with c ~ 1) in 
Euclidean space and time directions (though other options are possible). 

The matter and gauge fields live in these grid points and in the links that connect 
neighboring grid points, respectively. Either type of fields may have internal degrees 
of freedom (as indicated by spinor and/or color indices). For convenience one 
usually introduces periodic or antiperiodic boundary conditions. Hence, starting in 
{aNs, y, z, t) and hopping one unit in the 1-direction one ends up in (a 1, y, z, t) and 
similarly for the other directions, except for the 4-direction, where one picks up a 
sign if one is a fermion. With this setup the total number of degrees of freedom 
is finite, while the degrees of freedom themselves (the quark and gluon fields) are 
continuous (if they are discrete, too, one has a spin system). 

By introducing the rule that each configuration of quark and gluon fields shall 
be attributed the Boltzmann weight exp(— 5qcd), where the action |4] 

Nf 

Sqcd = ^Tt[F^uF^,) + g-(')(^+m«)g« + ie^e,,p,Tr{F,,Fp,) (1) 

^ i=l 

is the Euclidean counterpart of the Minkowskian Lagrangian, one defines the parti- 
tion function of QCD, provided there is an unambiguous measure for the variables 
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in each grid point or link [which is true for compact gauge groups Uke SU (Nc)]. This 
weight depends on the discretization or action chosen, i.e. on the details of how one 
expresses the quantities in ([I]) in terms of the quark fields q^^\x) and gauge links 
Ufj,{x), where x now denotes a 4-dimensional coordinate [which is implicit in 
For instance, the Dirac operator ]p may be anti-hermitian (as one would expect 
from a continuum viewpoint) in one discretization, but not even a normal operator 
in another one. Phenomenologically, the parameter 9 is known to be extremely 
close to (or tt), this is why the last term in ([T]) is usually omitted. 

To cut a long story short, the discretization of a field theory is an intermediate 
step which renders the path-integral (or partition function) finite. This happens 
right from the outset, in sharp contrast to perturbative approaches where one tames, 
at a much later stage, integrations which otherwise would produce infinities. The 
lattice provides a UV cut-off (through a > 0) and an IR-cut-off (through V = 
L^T < oo, where L = NgU is the box-length and T — Nta the time-extent) that 
respects the gauge invariance of Q. One then computes inverse correlation lengths 
in lattice units, e.g. a/^^^ = aM^^ and a/^o = gMq. The rule is that one must form 
ratios of such quantities before one is allowed to take the continuum limit. Hence, 
wile £,Tr/a would diverge, the ratio stays finite under a — 0. The crucial point 

is that the result is independent of the details of the lattice action ([!]), provided 
some guidelines were observed. In other words, the various lattice regularizations 
provide different ways of defining QCD, but they share a universal continuum limit. 
Any valid discretization of (IT]) is thus part of the definition of QCD ! 



2.2 Scale hierarchies in Lattice QCD 

With this bit of lattice ideology in place, it is time to reflect on something relevant 
in every days use of Lattice QCD, namely the scale hierarchies involved. Evidently, 
the shortest length-scale at hand is the lattice spacing a, the largest one is the box 
size L. The correlation lengths of the quarks that occur in ([T]) and of the hadrons 
that emerge as bound states should somehow lie in between these extremes. 

It turns out that for the shortest scales it is the (heaviest) quark mass that 
matters, while for the largest scales it is the mass of the (lightest) asymptotic state 
that is relevant. In other words, for the quarks we need to pay attention that they 
do not "fall through the grid" (i.e. we request aniq ^ 1), while for the pions we need 
to make sure that they are not squeezed too much (i.e. we request Mt^L 3> 1). This 
may be a bit of a challenge, because maintaining both conditions simultaneously 
requires having a large number of grid points in each direction. 

At the time of writing the largest number of grid points that can be sustained in 
current state-of-the-art lattice QCD simulations is of the order 128^. Even with this 
somewhat optimistic assumption on L/a the goal Mt^L > 4 would translate, at the 
physical pion mass, into L = 4.1/(135 MeV) = 6fm and hence into a = 0.047fm. 
With nic ~ 1.1 GeV the product is then aruc — 0.26. This number is small enough, 
but in order to take a continuum limit we need to have several lattice spacings. We 
can only afford to go to larger a; upon doubling the lattice spacing we end up with 
aruc — 0.52 which (with the actions currently in use) is barely acceptable. 

Hence simulating almost physical light quarks and, at the same time, charm 
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quarks without terrific cut-ofF effects poses enormous computationai requirements. 
Usually one chooses for the up and down quarks a common mass rriud^ and either 
this mass is kept in the vicinity of its physical (isospin averaged) value or nic is kept 
near its physical mass. Often some extrapolation in m„d and/or rric is performed, 
while nis is usually simulated next to its physical mass value. Simulating bottom 
quarks in the vicinity of their physical mass is out of question (except for special 
setups); they are usually treated in entirely different frameworks. And the top is 
omitted from Lattice QCD calculations, since it is too short lived to hadronize. 

2.3 How to remove systematic effects 

As emphasized in subsec. 2.1, the lattice (or something equivalent) is a necessary 
intermediate step in the definition of QCD. In consequence, a computer code which 
yields a stochastic estimate [S] of the path integral ([T]) for a given lattice spacing a 
and box size L, at a given set of quark masses niq = {rriudT'nis) [and possibly md, 
does not complete the job. A good lattice calculation combines the information 
from several simulations to remove the remnants of the lattice formulation. Thus, 
when reading a lattice paper one should have the following questions in mind: 

(1) Has the continuum limit (a— >-0) been taken ? 

(2) Are the finite- volume effects (from L < oo) under control ? 

(3) Are the simulations performed anywhere close to = 135 MeV ? 

(4) Advanced: are theoretical uncertainties properly assessed/propagated ? 

(5) Experts: algorithm details, treatment of isospin breakings, resonances, ... 

Unfortunately, all the interesting directions tend to be expensive in terms of com- 
puter time. The computational requirements tend to increase roughly like 

CPUcxl/a^-\ CPUcxL^ CPU cx l/mj^^ (2) 

and the meaning is that in each case the other parameters are held fixed. In reality, 
when pushing rriud down towards ~ 3.5 MeV also the box length L needs to 

be increased to maintain the bound AI-^L > 4. In addition, there are algorithmic 
issues which lead to a proliferation of noise near the chiral limit [6j; this makes 
simulations close to the physical mass point even more demanding. 

The quark fields in ([T]) may be integrated out, and this leads to a contribution 
^ log(det(^-|-m^*^)) to the effective action, where the sum runs over the flavors, 
i = 1, Nf. Todays lattice terminology attributes the label Nf = 2 to simulations 
with a common mass rriud of the "sea" quarks (i.e. those which come from the 
functional determinant). Studies which include a dynamical strange and possibly a 
dynamical charm quark are referred to as = 2-|-l and Nf — 2 -|- 1-1-1 simulations, 
respectively, to indicate that these quarks have separate masses. 

The lattice spacing and the quark masses cannot be dialed "a priori" because of 
renormalization effects. The simulations are governed by the bare gauge coupling 
/3 = 2Nc/g'^ and several (in the cases mentioned above: 2 or 2 -I- 1 or 2 + 1 + 1) 
bare mass parameters; the lattice spacing a and the quark masses niq are emergent 
quantities to be determined "a posteriori" . 
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3 Hadron Spectroscopy 

3.1 Quenched versus unquenched QCD 

Spectroscopy of stable states with a conserved quantum number (e.g. isospin) is 
about the easiest thing to do in Lattice QCD. One considers two-point functions 

Nf 

C{t) = (O(t)O^(O)) DU d{t)d\Q) TTdet(^[[/]-|-m(*') e"^!^! (3) 
Z J ^ ^ ' ^ ' 

"valence" "sea" 

where 0(.) is designed to absorb these quantum numbers, and at least one of the 
two involves some projection to a definite spatial momentum (e.g. p = 0). From the 
asymptotic behavior C{t) cx e~^'P-'* one extracts the mass of the particle. In the 
past the functional determinant was omitted {Nf = 0), while O would still contain 
so-called valence quarks; this is the quenched approximation. If a quark mass in the 
determinant is different from the mass of the same flavor in the interpolating fields 
0(.), i.e. TO^'Jf 7^ "^Md'' '^^^ ^^y^ ihaX the respective flavor is partially quenched. 



3. 2 Landscape of current Nf — 2 + 1 simulations 

As discussed before, generating ensembles with light vriud tends to be expensive. 
Fig.jl] depicts the sea pion masses and box sizes at which various collaborations 
managed to simulate. A significant number of ensembles is needed to extrapolate 
to (mP[j^',TOP'^y") and to remove the cut-off and finite-size effects. To reach the 
physical point most collaborations fix ~ m'^^^'^ and reduce ruud as much as 
possible (as of this writing only two collaborations can bracket m^J^ by the m^cL 
in the sea [8|9| ). By contrast, QCDSF pursues an interesting alternative. They 
start with a common sea quark mass near (2inn^^ -\-m^^^^) /Z and split the masses 
symmetrically such that this weighted sum stays invariant [10). 
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Figure 1. Landscape of simulated box sizes and pion masses in recent studies of full QCD. The 
amber shaded areas indicate the expected magnitude of relative finite size effects. Figure from [t]. 
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3.3 Reality check: spectra of stable hadrons 

We tend to postulate that the QCD Lagrangian ([T]) is the complete theory of strong 
interactions, valid both in the short-distance regime (where asymptotic freedom 
prevails) and in the long-distance regime (where confinement dominates). 

This is a highly non-trivial statement, and it was a historical achievement of the 
CP-PACS collaboration to show that quenched QCD (as an alternative candidate) 
does not pass the test [IT] (left panel in Fig. [2]). The unquenched version fares 
much better [l2] (TV/ = 2-1-1, right panel in Fig!j2]); of course, this does not amount 
to a proof. In the latter case three quantities are plotted with open circles, as they 
have been used to adjust m„(i, nig and to determine the individual a (cf. Sec. pi). 



3.4 Resonances and excited baryon spectra 

Those who wish to tackle harder problems may study how resonances interact with 
stable particles (note that the lattice results in Fig. [2] concern only their masses; 
the gray bands indicate experimental widths). A prominent example is gp 
coupling of the rho to two pions, on which there are nice results; see 13 
review. In the same category is the study of excited states with more mundane 
quantum numbers (e.g. those of the proton), see e.g. 14 for recent progress. 



TTTT ; the 

for a 



3. 5 Mixing of rj ~ rj' 

One of the long standing problems of QCD is to convincingly show that the mass 
splitting between the ry and rj' is due to the global axial anomaly (as argued by 
Witten and Veneziano long ago). This is now achieved on a quantitative level (i.e. 
with mixing) in two recent papers, one by RBC/UKQCD [iSj, one by ETMC [16] . 

3.6 Progress on glueballs 

Glueballs are hard in pure Yang-Mills theory (because they are noisy, and sometimes 
a vacuum contribution needs to be subtracted) , and they are extremely hard in the 
presence of dynamical flavors (since they mix with other flavor singlet states which 
do have valence quarks). For a recent paper on the subject see [Tt). 
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4 Hadron Structure 



4-1 Nucleon sigma terms and dark matter searches 

Dark matter experiments try to probe the interaction of unknown objects (WIMPs) 
with nuclei. This would happen via coupling to virtual light and strange quark loops 
in the nucleon, and the lattice is thus called to calculate the quantities 

o-^N o-ud = mud{N\uu+dd\N) , = ms{N\ss\N) (4) 

from first principles. This can be done either from the slope of Aljy versus and 
2M|- — (see Fig.|3]) or from a direct determination of the matrix elements. A 
straight average of all central values and total errors in Tab. [T] would suggest that 

CT^N - <yud = 47(9) MeV , a, - 48(25) MeV (5) 

is a conservative estimate (in some cases <Ts is obtained from un = "^mi/ms ■ Osjoud 
or /t, = m,{'N\ss\N) /Mn), see also [29; . With a-^N ^ d{MN)/d{M^) it foUows 
that the nucleon mass in the 2-flavor chiral limit is M^r (?Ti„rf = 0) ~ 880(20) MeV. 



4-2 Nuclear structure and recent news on qa 



The form factors of the nucleon with external 5", P, V, A currents and the respective 
radii are known to be hard on the lattice. Even the values at p^ = 0, in case of the 
axial current known as gA, turn out to be unexpectedly hard 30|31| (see also [32]). 
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Figure 3. The nucleon mass as a function of the average light quark mass ni u^ o c at fixed 
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4-3 Scattering ofm:, ttK, KK, ttN, NN and more 



Elastic scattering parameters can be extracted from the L-dependence of accu- 
rately measured energies of two-body states in a finite volume 33 , cf. Fig. |4] left. 
Recently a new method has been proposed to extract such information from two- 
body potentials, and it is claimed that such results are less sensitive to excited states 
contaminations ^ij. The right panel shows the central NN potential (spin-singlet 
channel) determined with this method; both the repulsion at short distance and 
the attraction at intermediate distance seem to grow as the pion mass decreases. 



4-4 From quarks to nuclei 

Evidently, there is a long way to go until all nuclear physics can be derived from 
Lattice QCD; see Fig.js] and 35|36 . A method to deal with the vastly growing 
number of contractions is proposed in 37 . The subject is reviewed in 38 . 
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in a finite volume (left, "level method" [33] , figure from [t]) and result for central nucleon-nucleon 
potential in ^So state at various pion masses (right, "potential method", figure from [34]). 
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5 Flavor Physics and FLAG 

5. 1 Light quark masses and ttstrong 

The quark masses and astrong are the parameters of the Standard Model that 
Lattice QCD can determine without further theoretical input. The starting point 
is a spectroscopy calculation like the one shown in the right panel of Fig. [2]. With 
this one has the quark masses in hand, albeit in an awkward scheme which is 
specific to the lattice action used in the course of the simulation. Before taking the 
continuum limit a conversion to a continuum scheme (e.g. MS or RI/MOM or SF) 
must be performed. A few recent computations are listed in Tab. [2]. 

Similarly, by observing how the lattice spacing a varies as one shifts the coupling 
/3 = 2Nc/g'^ or by matching small Wilson loops to lattice perturbation theory one 
can determine the QCD /3-function and the integration constant Aqcd (the case 
Nf =2 has been presented at this conference 44 ) or astrong (m) at a scale fj.^a~^. 



The PDG average of astrong (A^z) is dominated by the HPQCD results 45|40 . 



5. 2 Light decay constants and form factors 

The leptonic decay of a pseudoscalar meson [P — n, K, ...) is mimicked on the lattice 
by the coupling to an external axial current and captured in the decay constant fp. 
Experiment determines the width F c>c |/,rKi(iP or F cx l/ifVusP, and by dividing 
out fp one gets access to the CKM matrix element Vud or Vus, respectively. 

Similarly, the lattice can determine the — >■ tt transition form factor (with a 
vector current in between) f^'^'^{q^), and by combining it with the width of the 
semileptonic if — > tt decay one has another option to access V^s ■ 

In either case, when computing \Vud\'^ + iKisP + iKibP (the last contribution is 
tiny) one finds that the first-row CKM unitarity relation is well satisfied. 



5.3 Kaon mixing in and beyond the standard model 

The K^-K'^ oscillation is the main source of indirect CP violation. In the standard 
model (SM) it is dominated by two box diagrams which, via an effective field theory 
approach, may be shrunk into a ZXS" = 2 operator, such that one ends up with 
{K'^\Ovv+aa\K^) ■ The lattice can calculate this matrix element (which enters the 
unitarity triangle analysis) or Bk^ see e.g. |46|47|42]43 , as well as siblings which 



are relevant to beyond standard model (BSM) theories 48|49| 



Recently RBC/UKQCD managed to calculate the K — >■ (27r)/=2 amplitude 
(both Re and Im of this = 3/2 process) 50 1. Still, the M = 1/2 counterpart 
K — (27r)/=o (relevant to the M — 1/2 rule and e' /e) is significantly harder. 
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Table 2. Selection of recent computations of m^, m„d, "lu, "id [MeV] in the (MS, 2 GeV) scheme. 
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5.4 FLAG compilation and lattice averages 

The quantities discussed in the previous subsections have been evaluated by many 
lattice groups. In such a situation it is useful to have reviews which provide averages, 
like the one by FLAG [5l] (with a focus on low-energy physics) or [52] (with an eye 
on CKM physics). In either work the user is urged to cite the original literature. 
Recently these two teams merged into FLAG-II with the goal to provide a regularly 
updated compilation of a larger number of observables, with trustworthy assessment 
of the overall systematic uncertainty: http://itpwiki .unibe . ch/f lag| . 

5.5 Charm physics on the lattice 

Treating the charm quark on the lattice just became doable (cf. Sec.[2|. Recent years 
brought several precise determinations of (leptonic) decay constants and (semi- 
leptonic) transition form factors of D and Ds mesons [essentially the same story as 
was told for pions and kaons before]; see Fig. Island Tab. Is] for a tiny selection. 



5.6 Bottom physics on the lattice 



Due to its large mass the bottom quark tends to be treated with specialized methods 
(cf. Sec. [2]). Mass splittings in the T system and again decay constants and form 
factors are of interest; see Tab. |4] for a spotlight and 61|62 for recent reviews. 
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Figure 6. Shape of the form factors f^~'^(q'^) (normahzed at = 0.15 GeV^, left) and 
/_P~*^(g^) (normalized at = 0.10 GeV^, right) compared to experiment. Figure from Issl. 
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Table 3. Selection of recent computations of fjj [MeV], fr)^ [MeV] and fos/fo- 
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Recent Progress in Lattice QCD 



11 



6 Other Topics 

6.1 Large Nc, larger Nf, different representations 

Having a computer code which generates full QCD ensembles, it is straightforward 
to change the number of colors, Nc, the number of dynamical fermions, Nf, or even 
their color representation. Most practitioners are interested in the regime where 
such theories are slowly walking, i.e. near conformality, see 63|64 for reviews. 



6.2 QCD thermodynamics at fi = and fi > 

At vanishing baryon number density (/i = 0) and physical quark masses QCD shows 
a crossover [65|66| . Recently, a consensus was reached on the pseudocritical temper- 
atures Tc for a set of common observables (subtracted chiral condensate, Polyakov 
loop susceptibility) at physical quark masses and in the continuum 66|67|68 



QCD thermodynamics at non-zero chemical potential {fi > 0) is in a much less 
satisfactory state. The main reason is that there is a genuine sign problem in the 
regime >• T, Af^r (without any solution in sight). There are established techniques 
to investigate the regime /i ^ T, M^r (which is good enough for cosmology, since 
the strong transition was at /i ~ 0). To date it is not clear whether there is a line of 
first-order transitions separating the confined (hadronic) phase from the deconfined 
(plasma) phase. If it exists, it would terminate in a second-order endpoint [69|. 



6.3 Hadronic contributions to g—2 of the muon 

Hadronic contributions to vacuum polarization provide one of the major sources of 
systematic uncertainty in the computation of the anomalous magnetic moment of 
the muon = {g^ — 2)/2. The lattice can help by providing the (Fourier trans- 
formed) 2-point function of the electromagnetic vector current in such a hadronic 
environment. Recent papers include |70|71|72| ; see also 173]. 



6.^ QCD with isospin splitting and/or electromagnetism 

In current state-of-the art studies of lattice QCD photon loops are excluded from the 
simulations, and strong isospin breaking (due to m„ ^ mj) is ignored (since it is of 
the same order of magnitude). To compensate, input quantities (e.g. Mt^, Mk , Ms 
in Fig. [2] right) are usually adjusted for these effects, see [5l] for details. 

Recently, there has been significant progress in exploratory studies which include 
quenched photons (i.e. interacting with the valence but not with the sea quarks). 
The key issue is that the finite- volume effects are no longer exponentially small for 
large M^L, but only polynomially suppressed (see e.g. (74|). Similarly, there has 
been tremendous progress in splitting m„ ^ rud through reweighting techniques (see 
|75| , though in this work a considerable fraction of the reweighting power goes into 
shifting the average light quark mass downwards). 

The ultimate dream of todays lattice physicists is to perform Nf — 1 + 1 + 1 + 1 
simulations, such that the physical values of m„, m^, Wg, rric are bracketed by 
those in the ensembles, and the lattices would include the photon as a dynamical 
degree of freedom and cover, furthermore, several lattice spacings a and box sizes 
L such that an extrapolation a — and L — > 00 can be performed. 
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6.5 Towards exaflop machines 

Among numerically oriented scientists of all disciplines the lattice people are known 
for their hunger for computer time. Current flagship machines tend to provide 0(10) 
petaflops of peak performance [about 10^® multiply-and-add cycles per second] , with 
a power consumption of the order of 4MW. Reaching the exaflop era (10^® flops) 
requires new thoughts on how to make these machines even more energy efficient, 
since scaling the electricity bill up by two orders of magnitude is not an option. 

The linpack tests that are used to rank such machines on the "top 500" list use 
about 80% of the available cycles. For typical lattice codes this sustained perfor- 
mance ratio ranges between 20% and 50% (which is significantly higher than what 
is reached in most other disciplines). Still, with every new generation of supercom- 
puters the amount of parallelism grows, and it is an ever increasing challenge to 
adjust the code such that these high sustained performance figures would persist. 



6.6 More theoretical issues 

Let mc just mention some of the topics on which cither tremendous progress has 
been achieved or which mark open questions: improved actions and matching with 
perturbation theory, chiral symmetry in vector-like gauge theories, chiral gauge 
theories and CP violation, chiral symmetry and chemical potential, the generic sign 
problem at non-zero chemical potential, how one would formulate supersymmetry 
on the lattice, further particulars of the staggered fourth-root procedure, minimally 
doubled fermions and staggered fermions with non-standard mass terms, the issue 
of large autocorrelation times near the chiral and/or continuum limit. 



7 Summary 

Instead of giving a long summary, I try to formulate seven short messages: 

1. The discretization of the QCD Lagrangian ([ij [or something equivalent] is a 
necessary intermediate step in the definition of QCD. 

2. Spectroscopy of stable hadrons with 7V/ = 2orA^ = 2 + lorA^/ = 2 + l-|-l 
active (dynamical) quarks is a mature field. 

3. Spectroscopy of mixing or unstable [under strong interactions] states (often 
performed on the same gauge backgrounds) is developing fast. 

4. Lattice QCD yields vital input in CKM analysis and BSM bounds, much of 



which is reviewed in the compilations by FLAG 51 and latticeaverages 52 

Use these resources, but please be sure to cite the original papers ! 

There is rapid progress on a number of nuclear issues, such as the nucleon 

sigma term and various scattering lengths and phase-shifts. 

There is rapid progress on QCD thermodynamics. At /i = a consensus on 

the relevant pseudocritical temperatures Tc has been reached, while at /i 7^ 

it is still not clear whether a second-order endpoint exists. 

The lattice remains a relevant tool for conceptual work in field theory. For 
instance, it is still not clear whether supersymmetry is a well-defined field- 
theoretic concept (beyond perturbation theory). 
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